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Abstract: Wire-grid polarisers are versatile and scalable components 

which can be engineered to achieve small sizes and extremely high 
extinction ratios. Yet the measured performances are always significantly 
below the predicted values obtained from numerical simulations. Here 
we report on a detailed comparison between theoretical and experimental 
performances. We show that the discrepancy can be explained by the true 
shape of the plasmonic structures. Taking into account the fabrication 
details, a new optimisation model enables us to achieve excellent agreement 
with the observed response and to re-optimise the grating parameters to 
ensure experimental extinction ratios well above 1,000 at 850 nm. 
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1. Introduction 

Although the concept of microwave polariser based on subwavelength gratings dates back to 
the 19^^ century, it could not be adapted to the optical domain until the development of mi¬ 
crofabrication technologies in the 1960s Q. Since then, these so-called Wire-Grid Polarisers 
(WGP) have been used frequently as they offer small footprints and large acceptance angles, 
and can now be found in daily-life devices such as sensors O or displays. Current research 
focuses onto low-cost and high volume manufacturing processes, UV applications (O, as well 
as high performance. 

A WGP consists of parallel metal stripes on a transparent substrate, such that the period p 
is smaller than the wavelength of the incident radiation. In this configuration, the Transverse- 
Magnetic wave (TM, also 7t or p), polarised orthogonally to the stripes, can experience Extraor¬ 
dinary Optical Transmission (EOT) due to coupling to surface plasmon polaritons (H O and 
waveguiding effect through the slits Ena. The subwavelength dimensions guarantee zero-th 
order diffraction. On the contrary, the Transverse-Electric polarisation (TE, also called cr or 
s) is almost perfectly refiected. This leads to high polarising efficiencies, characterised by the 
Extinction Ratio, i.e. the ratio of the transmission of orthogonal polarisations ER = Ttm/Tte- 

Usually a trade-off has to be found between high ER and high transmittance of the TM 
modes. Our application for quantum cryptography [0 requires an array of four equivalent po¬ 
larisers with a minimum ER of 1,000 at 850 nm, while the overall transmission is not a critical 
parameter. Numerical simulations typically used to find this compromise promise very high ER 
which are, however, never measured in real devices. Discrepancies of several orders of magni¬ 
tudes have been observed Enniini . To achieve better experimental results, new structures such 
as dual-gratings have been proposed CD. 

Here we show that geometrical deviations of the manufactured structure with respect to those 
used in the simulations are the reason for this discrepancy IT^ . In particular our gold grat¬ 
ings exhibit trapezoidal stripes instead of rectangular ones, a common problem, with only a 
few exceptions |[T3l[T4]|, to all fabrication techniques such as anisotropic etching, Eocused Ion 
Beam (EIB) milling and nanoimprint. We observe an approximately exponential decrease of 
the extinction ratio with the tilting angle of the slits. Yet accounting for the real shape of the 
wires, in particular for the average slit width, leads to very good agreement with the observed 
performances. Our improved numerical simulation model predicts that high extinction ratios 
can be achieved even with imperfect structures. By reducing the design slit width we obtain an 



experimental ER of better than 1,000 and a transmission of 9 % for a wavelength of 850 nm. 

This paper is organised as follows: we first compare the grating parameters and the simulation 
results based on perfectly rectangular structures (Section 2) with the performance of samples 
fabricated accordingly, exhibiting high discrepancies with the numerical computation (Section 
3). In Section 4 we describe a refinement of our simulation model and an optimisation step to 
achieve the desired performances. Finally we confirm our new model with the characterisation 
of new optimised samples. 

2. Standard grating optimisation 

Considering a gold grating on a glass substrate, the parameters to be optimised are the period 
p, slit width w and thickness /z, as depicted in Fig.[2 The period as well as the duty cycle (ratio 
w to p) are directly related to the coupling of the TM polarisation with surface plasmon on the 
top and on the bottom of the metal stripes. Certain combinations of parameters lead to efficient 
field confinement between the stripes (Rayleigh-Wood resonances), thereby increasing the 
resulting ER. In general smaller slits tend to a favourable increase of the polarisation selection 
for a fixed period, but also lead to lower TM transmission. The latter is nonetheless also depen¬ 
dent on the metal height, which should be selected carefully to guarantee efficient tunnelling 
through the slit via mode-matched Fabry-Perot cavity resonances. Increasing thickness leads 
however to an exponential decay of the transmission of the TE polarisation, mainly responsible 
for the evolution of the ER. The interplay between both vertical and horizontal resonance 
phenomena is studied in (TS) [Tbl. Yet due to fabrication limitations, high quality high as¬ 
pect ratio structures are hard to manufacture, imposing a constraint on the grating parameter set. 


TE TM 



Fig. 1: Geometry of a gold wire-grid polariser on a glass substrate. 


Previous studies El indicate typical grating parameters, and a period of ;? = 500nm seems 
to deliver good performances around A = 850 nm. In this regime the condition p <C A is not 
fulfilled, therefore the Effective Medium Theory (EMT) cannot be used. The optimisation of 
the slit width as well as the gold thickness was thus studied using FDTD simulations (MEEP 
ca). This simulation program uses the Lorentz-Drude dispersion model with optical constants 
obtained from thin films measurements El . Due to the very long computation time associated 
with a resonant structure, the spatial resolution was limited to 7 nm. Figure |2(a)| shows the 
dependence of the extinction ratio on the thickness and the slit width. Beyond h = 320 nm, the 
ER reaches very high values but we rather concentrate on the first vertical resonance at /z = 
270 nm, as the fabrication of very narrow slits with straight fianks becomes more difficult with 
increasing metal thickness. We also observe on Fig. |2(b)| a horizontal resonance at w = 90nm, 
where the ER increases due to higher transmission of the TM-modes. 
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Fig. 2: Influence of the design parameters on the performance of the polariser. (a) Dependence 
of the extinction ratio on h and w for /? = 500 nm. Here the effect of the thickness is considered 
only around the first vertical resonance. The green line represents our design rule of ER = 1000. 
(b) Transmission of the TM-modes as a function of w for /z = 270 nm. 


3. Fabrication and characterisation 

A 170 pm thick glass substrate was coated with 3nm Titanium (adhesion layer) followed 
by 265 nm gold using Electron Beam Physical Vapour Deposition. Deposition rates were 
maintained relatively low (below 1.5A • to minimise the granularity and the surface 
roughness of the layers. Although simulations indicate a resonance for h = 270nm, better 
performances were experimentally observed for h = 265 nm. This shift was also predicted by 
additional high-resolution simulations taking into account the presence of the adhesion layer. 
The theoretical ER obtained for the experimental configuration is nevertheless comparable with 
the one resulting from the simple model of gold Qi = 270 nm) on a glass plate. The sub wave¬ 
length gratings (120 x 120 pm) were then engraved using the Focused Ion Beam (FIB) milling 
technique (Zeiss Cross-Beam) at the Center for Nanostructures and Nanomaterials (ZNN, 
Munich, Germany). The wire-grid polarisers were characterised by measuring the ER for well- 
defined input polarisations. While this study focuses on performance at 850nm, we refer the 
reader to (m for similar numerical and experimental results including wavelength dependence. 


Table 1: Comparison between simulated and experimental extinction ratios for different slit 
widths. Here the structure is simulated with a period p = 500 nm a gold thickness of /z = 270 nm, 
but fabricated with h = 265 nm with additional 3 nm Ti. 


Sample w (nm) ER (experimental) ER (simulated) 


Ax 

120 

380 

Ai 

120 

650 

Bx 

80 

720 

Bi 

80 

850 


2520 


12700 


According to the simulations, an extinction ratio above 2,000 should be obtained for slit 
widths below 120nm. The first sample was fabricated with this largest possible width (Samples 
A, see Table [^, yet the measured values did not exceed 1:650. In order to achieve the desired 








































performances, new samples with smaller slit widths were manufactured. Although the experi¬ 
mental ER slightly increases for w = 80 nm, the simulations indicate a much stronger increase 
by at least a factor of 10. Moreover, the measured ER shows a large scatter for different samples 
although the gratings exhibit a similar geometry and roughness in top view SEM pictures, typ¬ 
ically similar to Eig. |3(a)| A significant difference appears when investigating the cross-section 
of the gold stripes using EIB milling. The side-view SEM pictures in Eig. |3(b) now reveal a 
rather trapezoidal profile, with a non negligible difference between the lower and upper base 
lengths of about lOOnm. This shape results from side redeposition of the ablated material during 
the milling of the stripes, which can be potentially reduced by using a multi-pass EIB writing 
mode. 



Eig. 3: (a) Top and (b) side-view SEM images of the gold stripes, (c) Geometrical model taking 
into account the trapezoidal shape of the wires. 


4. Grating optimisation with realistic geometrical model 

In order to understand the often reported reduced performances of the manufactured WGP 
compared to predictions, a deeper study involving further simulations has been carried out. 
To account for the observed geometry we introduce a refined model based on a more realistic 
trapezoidal stripe shape. The parameter w now denotes the slit width at the bottom of the stripes, 
close to the substrate, and the angle a corresponds to the deviation from a perfect rectangular 
structure, as illustrated in Eig. |3(c)| 

The transmission of both TE and TM modes, as well as the ER are computed varying again 
the slit width w and additionally the opening angle a. Eigure [4(^ shows an exponential de¬ 
crease of ER when varying the angle a, but also that the trapezoidal shape leads to slightly 


higher transmissions for the TM modes, see Eig. 4(b) This is consistent with the case for rect¬ 
angular shape, where Tjm increases with w. This effect is nonetheless reduced in the case of 
the resonance (w = 90nm). While the influence of a on TM-states is limited, it has dramatic 
consequences for the TE-polarisation, as can be seen in Eig. |4(c)| In this case we observe an 
exponential increase of the transmission with the opening angle, which consequently leads to 
the exponential decrease of the ER. When using the trapezoidal shape in the simulations, the 
theoretical ER values are reduced from 12700 (rectangular stripes, w = 120nm) down to 500 
for B\ {oCexp = 25^) and 2700 for B 2 {(Xexp = 16^). This trapezoidal shape indeed explains the 
order of magnitude discrepancy observed between simulations of perfectly rectangular stripes 
and realistic samples. 

Since Ttm is largely independent of the angle a, it is the dependence of the transmission of 
the TE-mode Tte which determines the ER. A closer inspection adopting the model of guided 
waves in the slit 1^ shows that the TE-mode is exponentially damped as the wavelength of 


the incoming field is above the cut-off wavelength. Eigure |5(a)| shows exemplarily the result 
of the numerical simulations for two different heights. Evidently, the field is exponentially 
damped when propagating through the slit, which also explains the exponential dependence of 
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Fig. 4: Dependence of the performance on the slit angle a (h = 270nm): (a) extinction ratio, 
(b),(c) transmission of TM and TE polarisation modes, respectively. The exponential increase of 
the transmission of TE-modes is mainly responsible for the significant reduction of the extinc¬ 
tion ratio, (d) Transmission of the TE polarisation as a function of the effective slit width w. 


the transmission on the metal thickness. Analysing the dependence of the damping coefficient 
on the slit width shows an approximately linear behaviour in the region of relevant slit widths 
in Fig. |5(b)[ which is also in good agreement with calculations of effective refraction indices 
in previous studies (T\} . Given smooth variation of the slit width, Tte and thus ER can be well 
approximated taking the effective, i.e., averaged, width of the slits into account. Figure [4(d^ 
shows that the unnormalised values of Tte for different average slit widths w = {wmax'^'^min)/^ 
are indeed approximately independent of the angle a. 

Using these findings for the next optimisation step, we note that the targeted extinction ratio 
of 1,000 can be achieved for rectangular stripes with w < 1 lOnm, but not for achievable angles 
a > 16^. Yet, reducing the slit width down to 70nm (corresponding to w = llOnm) brings a 
clear improvement. This simulation was verified experimentally by the fabrication of a new po- 
lariser array, presented in Fig.[^ Table [^presents the characterisation of these four samples and 
compares them to the theoretical performances simulated with the refined model. As expected, 
the ER exceeds 1,000 and even reaches 1,800, yielding the best ER observed so far for 850 nm. 
The transmission is similar for all samples and reaches 9%. The simple simulation model based 
on trapezoidal stripes is therefore suited for realistic devices as it exhibits excellent agreement 
with the experimental data. 

The remaining difference between experiment and theory is most likely due to the spatial 
resolution of the simulation and the measurement uncertainty of a on the SEM image, as well 





















































Fig. 5: (a) Simulated intensity distribution for rectangular stripes separated by w = 70nm and 
w = 130nm respectively. The small insets show the exponential damping of the intensity within 
the slit (I{z) ^ (b) Damping constant 7 extracted from the intensity profile for different 

slit widths. 


as the aforementionned effect of the adhesion layer. According to recent studies surface 
roughness should not degrade the performance significantly. In addition, the trapezoidal shape 
smoothens the field distribution compared to rectangular shapes, further reducing the influence 
of surface irregularities. 



Fig. 6 : SEM pictures of a four-polariser array exhibiting extinction ratios up to 1,800. (a) Top 
view of the matrix, (b) Cross-section of the fourth grating. The decrease in performances asso¬ 
ciated with a large tilting angle a > 16^ was compensated by reducing the slit width. 


Table 2: Experimental results obtained after optimisation of the geometry using a trapezoidal 
model. The data exhibit clearly improved agreement with the simulations. 


Sample 

w (nm) 

Wmax (nm) 

aO 

ER (experimental) 

ER (simulated) 

1 

70 

150 

16 

1800 

4200 

2 

70 

160 

19 

1620 

2870 

3 

80 

160 

16 

1200 

1510 

4 

70 

175 

21 

1150 

1544 


5. Conclusion 

We presented clear evidence of the impact of geometrical deviations from perfect rectangular 
cross-section onto the performance of wire-grid-polarisers. Transversal SEM pictures indicate a 
trapezoidal shape of the stripes, a defect present with all currently used fabrication techniques. 































While better rectangular stripe profiles could be achieved, there will always be some imper¬ 
fections left. Nevertheless, the real performance of wire-grid polarisers can be well simulated 
when accounting for the true shape of the wires or in a first approximation for the effective slit 
width. High extinction ratios can be achieved even with imperfect structures when optimising 
other parameters. Comparison between experiments and simulations shows, for the first time, 
significantly improved agreement, with experimental devices reaching extinction ratios of up 
to 1,800 and transmission of 9% at 850nm. 
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